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THE  REASON  FOR  PERFORMING  THIS  WORK  was  to  improve  the  accuracy  and  efficiency  of 
the  target  acquisition  representation  in  the  Nuclear  Fire  Planning  and  Assessment  Model,  Version  III 
(NUFAM  in),  and  develop  a  methodology  for  estimating  the  sensitivity  of  the  model  output  to 
inputs  related  to  acquisition  and  movement. 

THE  PRINCIPAL  FIN  DINGS  of  the  work  reported  in  this  paper  follow. 

(1)  The  target  acquisition  and  unit  movement  processes  may  be  represented  using  alternating 
renew^  processes. 

(2)  The  current  data  for  target  acquisition  used  within  NUFAM  III,  the  probability  of  target 
acquisition  (POTA)  values  for  a  2-hour  glimpse  period  derived  from  the  Target  Acquisition  Study, 
Version  HI  (TAS  III),  require  the  number  of  periods  until  acquisition  to  be  distributed  as  geometric. 

(3)  The  process  of  tracking  nuclear  targets  after  acquisition  may  be  represented  based  on  POTA 
values  if  certain  assumptions  are  accepted. 

THE  MAIN  ASSUMPTIONS  follow. 

(1)  The  POTA  values  derived  from  TAS  HI  may  be  used  to  represent  the  probability  of 
acquisition  during  successive,  independent  2-hour  glimpse  periods. 

(2)  The  uncertainty  about  the  probability  of  continuing  to  detect  during  a  2-bour  glimpse 
period,  given  that  the  target  has  been  previously  detected  (i.e.,  the  probability  of  tracking)  can  be 
described  using  a  uniform  (  POTA,  1  )  distribution. 

(3)  Units  are  acquired  and  move  independently  of  each  other. 

THE  PRINCIPAL  LIMITATION  of  this  paper  is  the  reliance  on  the  POTA  values,  which  do  not 
provide  an  accurate  representation  of  the  time-dependent  probability  that  a  unit  is  on  the 
acquisition  list  at  any  given  period  of  time.  Nevertheless,  the  POTAs  represent  the  best  available 
data  until  new  acquisition  and  tracking  distributions  are  developed  from  the  Target  Acquisition 
Study  rv  (TAS  FV)  study. 

THE  SCOPE  OF  THE  TECHNICAL  PAPER  is  limited  to  representing  the  target  acquisition 
process  in  the  March  1989  version  of  NUFAM  HI  using  the  current  input  data  values  and  formats. 

THE  TECHNICAL  OBJECTIVES  were  to  compute  the  probabilities  that  a  unit  in  NUFAM  HI: 

(1)  is  acquired  and  remuns  on  the  acquisition  list, 

(2)  is  available  for  fire, 

(3)  has  moved  at  the  time  of  detonation, 

(4)  can  be  fired  at  by  a  nuclear  weapon,  and 

(5)  is  hit  by  a  nuclear  weapon. 

THE  BASIC  APPROACH  used  alternating  renewal  processes  to  represent  the  acquisition  status  and 
movement  state  of  each  unit  within  NUFAM  III. 


t; 


THE  TECHNICAL  PAPER  SPONSOR  was  the  US  Army  Concepts  Analysis  Agency. 

THE  TECHNICAL  PAPER  was  written  by  MAJ  Mark  A.  Youngren,  Requirements  Directorate,  US 
Army  Concepts  Analysis  Agency. 

COMMENTS  AND  QUESTIONS  may  be  addressed  to  the  Director,  US  Army  Concepts  Analysis 
Agency,  ATTN:  CSCA-RQN,  8120  Woodmont  Avenue,  Bethesda,  MD  20814>2797. 
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CHAPTER  1 
BACKGROUND 


1-1.  INTRODUCTION.  The  US  Army  Concepts  Analysis  Agency’s  (CAA)  Nuclear  Fire  Planning 
and  Assessment  Model  III  (NUFAM  III)  explicitly  simulates  the  acquisition  of  mobile  nuclear  targets 
in  order  to  generate  an  acquired  unit  set  for  targeting  at  the  time  of  fire  planning.  The  version  as  of 
March  1989  relies  on  input  values  (defined  below)  such  as  the  POT  A,  retention  time,  travel  time, 
and  stay  time.  This  paper  addresses  an  analytic  representation  of  the  target  acquisition  model  ustng 
the  data  and  assumptions  contained  in  the  current  version  of  NUFAM  III.  As  a  result,  it  does  not 
incorporate  more  advanced  methods  for  modeling  the  stay,  move,  acquisition,  and  retention  times.  It 
is  intended  to  be  used  in  support  of  NUFAM  III  —  based  analyses  (such  as  the  NUX-97  study)  and 
in  support  of  the  nuclear  effects  module  (NEMESIS)  of  the  Integrated  Warfare  Force  Evaluation 
Model  (IWFORCEM),  until  improved  data  and  methodology  are  available  from  the  TAS  IV  (Target 
Acquisition  Study  IV)  and  TAME  (Target  Acquisition  Methodology  Enhancement)  Study  (Bauman 
[1989]).  A  basic  knowledge  of  the  terms  associated  with  the  NUFAM  III  model  and  input  data  base 
is  assumed.  A  complete  explanation  may  be  found  in  the  NUFAM  III  User’s  Manual  (Schuetze  and 
Albrecht  [1986]). 

1-2.  THE  TARGET  ACQUISITION  PROCESS 

a.  The  process  of  target  acquisition  is  defined  as  the  process  that  searches  for.  detects,  and 
identifies  potential  target  units  for  placement  on  an  acquisition  list.  A  tinif  is  a  military  organization 
composed  of  personnel  and  equipment  {target  elements)  that  is  identified  in  NUFAM  III  as  a 
potential  target.  All  units  are  categorized  into  type  units  according  to  their  characteristics  (types  of 
equipment,  movement  rates,  etc.).  There  are  usually  many  units  of  each  type  to  be  found  within  the 
NUFAM  III  data  base.  A  unit  is  acquired  when  it  is  placed  on  the  list;  acquisition  is  lost  when  it  is 
removed  from  the  list.  A  unit  may  be  dropped  (removed)  from  the  list  either  due  to  a  negative 
sensor  report  (i.e.,  we  no  longer  detect  its  presence  when  looking)  or  it  may  be  rlronned  after  some 
period  of  time  when  the  acquisition  information  cannot  be  updated.  The  time  to  acquisition,  T^,  is 
the  time  it  takes  to  acquire  a  unit  once  any  previous  ^u:quisition  has  been  dropped.  The  retention 
time,  T^,  is  the  length  of  time  the  unit  is  retained  on  the  acquisition  list. 

b.  Targets  are  units  that  are  planned  for  engagement  by  tactical  nuclear  weapons.  Relocatable 
targets  are  units  which  have  the  capability  to  move  during  the  scenario  of  interest  (although  they 
may  or  may  not  retain  mission  capability  during  movement).  As  a  consequence  of  this  movement 
capability,  units  do  not  remain  acquired  indefinitely  (unless  they  can  be  tracked  indefinitely  once 
acquired);  at  some  time,  they  move  and  the  acquisition  is  lost.  Even  if  a  monitoring  or  tracking 
capability  exists,  there  is  a  probability  that  such  tracking  will  be  lost  over  time.  NUFAM  III  has 
recently  been  modified  to  include  fixed  targets  as  well  as  relocatable  targets,  but  the  fixed  targets 
have  predefined  locations  and  are  considered  to  be  “acquired”  at  the  beginning  of  the  simulation 
and  retained  indefinitely  —  thus  the  various  measures  discussed  in  this  paper  do  not  apply. 

c.  Units  that  are  engaged  using  conventional  weapons  are  generally  fired  upon  soon  after 
acquisition.  Nuclear  targets,  however,  differ  from  conventional  as  they  are  planned  for  specific 
purposes  dictated  by  the  overall  theater  situation.  As  a  result,  nuclear  targets  are  not  normally 
engaged  as  they  are  acquired;  rather,  nuclear  fires  are  directed  at  targets  that  are  acquired  and 
perceived  to  still  be  in  place  at  the  time  that  the  weapons  are  approved  for  fire. 
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d.  Figure  1-1  illustrates  a  representative  acquisition  sequence  for  a  relocatable  unit.  Any  given 
unit  will  be  placed  in  only  one  of  two  states:  acquired  or  not  acquired.  Once  the  unit  has  been 
dropped  from  the  list,  it  is  immediately  subject  to  being  reacquired.  A  target  acquisition  process  is 
therefore  a  temporal  series  of  such  acquisition  states.  Since  nuclear  weapons  use  will  generally  occur 
only  after  the  conventional  battle  has  been  underway  for  some  time,  we  are  interested  in  the  target 
acquisition  probabilities  at  some  point  in  time  well  after  target  acquisition  has  begun.  For 
relocatable  targets,  the  target  acquisition  process  can  be  approximated  as  an  alternating  renewal 
process  of  indefinite  length.  Both  the  time  that  the  target  is  dropped  from  the  list  and  the  time  that 
the  target  is  acquired  ue  renewal  points  of  this  alternating  renewal  process. 


Acquired 
Not  Acquired 


Time 


■> 


Figure  1-1.  Possible  Target  Acquisition  Sequence 


1-3.  THE  FIRE  PLANNING  AND  EXECUTION  SEQUENCE  IN  NUFAM  III.  NUFAM  III 
represents  tactical  nuclear  warfare  as  a  series  of  engagements.  For  each  engagement,  each  side  has  a 
fire  planning  cycle  of  fixed  duration  during  which  nuclear  fires  are  planned  against  available  units  on 
the  acquisition  list.  The  fire  planning  cycle  terminates  at  the  beginning  of  some  specified  time 
interval  within  which  weapons  may  be  fired,  which  is  referred  to  as  the  fire  execuiton  cycle.  The 
time  that  the  fire  planning  cycle  ends  is  denoted  as  t^;  the  NUFAM  III  model  assumes  that  any  unit 
on  the  acquisition  list  at  time  tp  is  available  for  fire  planning.  A  unit  that  may  be  planned  as  a 
target  will  have  had  a  recent  acquisition  time  Tacf  prior  to  tp.  The  nuclear  round  will  detonate  at 
some  time  T,j  during  the  fire  execution  period  (Figure  1-2). 


Fire  Planning  Cycle 


Fire  Execution  Cvcle 


Time 


Figure  1-2.  Acquisition,  Planning  and  Detonation  Events 


1-4.  SOME  CONCEPTS  BEHIND  ALTERNATING  RENEWAL  PROCESSES 

a.  Representing  the  target  acquisition  process  as  an  alternating  renewal  process  allows  us  to  use 
the  extensive  theory  developed  around  such  processes.  If  we  enter  such  a  process  for  a  unit  at  a  point 
T  in  time  chosen  independently  of  the  target  acquisition  sequence  of  that  unit,  the  point  T  is  more 
likely  to  lie  in  a  long  interval  than  a  short  interval.  This  fact  needs  to  be  kept  in  mind  when 
determining  the  distribution  of  the  time  remaining  from  our  entry  point  T  until  the  next  transition 
(acquisition  or  drop  event),  or  the  time  since  the  last  transition  to  the  time  T.  Fortunately,  these 
distributions  have  been  worked  out  (cf.  Ross  [1983]).  Following  the  terminology  common  in 
reliability,  I  call  the  time  from  the  last  transition  to  our  random  entry  point  the  age  and  the  time 
from  our  random  entry  point  to  the  next  transition  the  residual  life.  An  e.xample  using  the 
acquisition  sequence  showing  the  age  and  residual  life  of  the  distribution  of  the  retention  time  T^^  is 
shown  in  figure  1-3. 
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Acquired 

Not  acquired 
Time 

T 

Figure  1-3.  Age  and  Residual  Life  of  the  Retention  Time  Distribution 

b.  Using  the  retention  time  T^  cis  an  e.xample,  the  random  variable  for  the  time  remaining  from 
our  entry  point  T  until  the  next  transition  (drop  event),  which  is  the  residual  life  of  the  retention 
time,  is  denoted  as  Y^.  The  general  form  for  a  distribution  of  a  residual  life  is: 

P(  Y«  >  j, )  =  1  -  |%(  T;j  >  u  )  du,  [1] 

where  E[T^]  denotes  the  mathematical  expectation  of  the  random  variable  T^  (the  expectation  may 
be  thought  of  as  the  “average”  of  T^).  Since  T^  is  fixed  a.  constant  value  R  in  NUFAM  III. 
E(T/j]  =  R  and 

F.J.  (u)  =  P[  T;j  >  u  ]  =  1,  0  <  u  <  R,  0  for  u  >  R.  Thus 

R 

P(  Y;j  >  y  )  =  1  -  ^  ^  1  du  =  I  -  I,  0  <  y  <  R;  0  for  y  >  R.  [2] 

1-5.  CURRENT  TARGET  ACQUISITION  MODELING  IN  NUFAM  III 

a.  The  target  acquisition  process  presently  modeled  in  NUFAM  III  relies  on  a  probability  of 
operational  target  acquisition  (POTA),  provided  by  the  Target  Detection  Routine  (TADER)  model 
as  part  of  the  TAS  III  study  (Penn  and  Bauman  [1987]),  as  a  measure  from  which  the  time  to 
acquisition,  T^,  can  be  computed.  Given  many  characteristics  of  the  sensors,  target  units,  and 
battlefield  scenario,  the  POTA  represents  “the  probability  of  detecting,  recognizing,  and  locating 
various  types  of  potential  targets  at  prescribed  distances  from  the  forward  line  of  own  troops 
(FLOT)  during  a  random  but  limited  period  of  time  in  a  day  of  intense  combat" 
(Bauman  and  Penn  [1987]  ),  where  the  period  of  time  is  currently  set  at  2  hours.  The  POTA  varies 
by  type  unit  and  zone,  where  the  zone  defines  a  FLOT-to-target  range  band;  Zone  1  is  0-3  km.  Zone 
2  is  3-12  km.  Zone  3  is  12-25  km.  Zone  4  is  25-100  km,  and  Zone  5  is  100-300  km. 

b.  In  order  to  apply  the  POTA  to  .x  time  period  larger  than  2  hours,  NUFAM  III  assumes  that 
the  probability  of  detecting  a  unit  during  any  .i-hour  period  is  independent  of  the  probability  of 
detecting  that  unit  during  any  different  2-hour  period.  Thus  the  probability  of  first  detecting  a  unit 
during  the  i-l-Ist  2-hour  time  p)eriod  is  equal  to  POTA  •  (  1  —  POTA  )  ,  /t  =  0,  I,  .  ..  .  which  is  a 
simple  geometric  distribution.  NUFAM  III  draws  against  the  POTA  every  2  hours,  which  is 
equivalent  to  drawing  from  this  geometric  distribution.  Once  NUFAM  III  has  determined  the  2-hour 
period  within  which  an  acquisition  is  to  occur,  it  schedules  the  actual  time  of  acquisition  randomly 
within  the  2-hour  interval.  Thus,  if  a  random  variable  X  is  defined  distributed  as 
geometric  (POTA),  and  a  random  variable  U  is  defined  distributed  as  uniform  (0.2).  then  the 
random  variable  for  the  time  to  acquisition.  T^,  is  identically  distributed  as  (thus  can  be 
represented  as)  2X  +  U. 
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c.  The  retention  time,  T^,  varies  with  the  capabilities  of  the  sensor  systems  to  revisit  or  track 
previously  acquired  units.  Various  sensor  system  sceDar;'>s  are  discussed  in  Chapter  2.  The  retention 
time  is  currently  modeled  in  NUFAM  III  os  a  fixed  input  value  R  which  represents  the  expected 
value  of  T^;  at  some  future  time,  it  may  be  desirable  to  model  the  retention  time  as  a  random 
variable  for  some  sensor  system  scenarios. 
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CHAPTER  2 

ANALYTIC  REPRESENTATION  OF  TARGET  ACQUISITION 


Section  I.  INTRODUCTION 


2-1.  RENEWAL  REPRESENTATION  OF  THE  TARGET  ACQUISITION  PROCESS 

a.  It  is  not  necessary  to  explicitly  represent  the  target  acquisition  process  in  NUFAM  III.  In 
NUFAM  III,  we  are  interested  in  two  things  related  to  the  acquisition  process:  the  probability  that  a 
unit  is  on  the  acquired  unit  list  at  the  time  that  we  conduct  our  nuclear  fire  planning,  and  the 
probability  that  the  unit  will  remain  on  the  acquisition  list  until  engagement  by  a  nuclear  weapon. 
Representing  the  target  acquisition  process  as  an  alternating  renewal  process  allows  us  to  easily  find 
the  probability  of  acquisition  at  any  random  point  of  time,  pae^: 

P"*’  “  ^  Youngren  [1988b]  ),  [3] 

b.  In  NUFAM  III,  T^  =  2X  -H  U  and  E[T^]  =  R,  where  R  is  a  fixed  retention  time  for  th> 
type  unit  and  zone  (recall  that  X  —  geometric  (POTA)  and  U  ~  uniform  (0,2)  ).  Thus  for  any 
given  unit, 

R  R 

P*'*  =  2E[X]  -b  E[U1  4-  R  =  „(l-POTA)  .  ,  .  „ 

2  ■  POTA  ”  +  1  R 


~  — 5 - ; — ^ 

pQta  ~  ^ 

c.  Example.  Suppose  that  the  POTA  IS  0.20  and  the  retention  time  R  is  2  hours.  Then  E[  T^  ] 
=  El  2X  1  +  El  U  1  =  2  1  =  2  ("f  +  1  =  9  hour.  »d 

~  2-4  +  I  +  2  ~  n  ~ 


2-2.  RELATIONSHIP  TO  THE  STAY-MOVE  CYCLE 

a.  NUFAM  III  simulates  movement  by  causing  the  unit  to  vary  between  two  alternating  states: 
the  move  state  (moving),  or  the  stay  state  (stationary).  We  define  the  variable  S  to  represent  the 
average  length  of  time  that  a  unit  is  stationary,  each  time  that  it  is  stationary  in  the  scenario,  and 
the  variable  M  to  indicate  the  average  length  of  time  that  it  is  moving,  each  time  that  it  is  moving 
in  the  scenario.  These  values  S  and  M  are  the  best  estimates  of  the  average  stay  and  move  times  for 
each  type  unit  represented  in  the  model  scenario.  The  random  variables  representing  the  stay  and 
move  times  generated  by  the  model  are  denoted  as  S^f  and  Mj\f,  respectively.  The  unit  movement  is 
represented  as  an  alternating  renewal  process  just  like  the  target  acquisition  process  (Figure  2-1). 
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Figure  2-1.  Possible  Target  Acquisition  and  Stay/Move  Sequences 


b.  The  values  input  into  the  NUFAM  III  Type  Unit  file  are  the  “average  stay  time”  denoted  as 
S^  and  the  “travel  time”  denoted  as  for  each  type  unit.  In  some  cases,  we  need  to  provide 
NUFAM  III  with  numbers  that  do  not  reflect  the  scenario  in  order  for  the  probability  that  a  unit 
will  be  hit  by  a  nuclear  weapon  to  be  computed  correctly  within  the  model.  In  these  cases,  we 
distinguish  between  the  average  stay  and  move  time  numbers  input  to  NUFAM  III,  Sj^^  and  Mj^^, 
with  the  scenario  driven  average  stay  and  move  times  S  and  M. 


c.  A  target  can  be  acquired  and  observed  while  either  stationary  or  moving,  so  T^  and  T^  may 
be  dependent  upon  the  distributions  of  the  scenario  stay  and  move  times  S  and  M.  However,  in 
NUFAM  III  we  assume  that  if  the  parameUrs  of  the  distributions  of  T^,  T^,  S  and  M  are  known, 
then  the  variables  T^,  T^,  S  and  M  are  mutually  independent.  Thus  we  can  examine  the  unit 
move/stay  state  at  any  arbitrary  point  in  time  independently  of  the  target  acquisition  state.  Within 
the  NUFAM  III  model,  the  probability  that  a  unit  is  in  the  stay  state  (stationary)  at  any  arbitrary 
point  in  time  is  p,<4,  : 


Pjtajr 


[5] 


where  is  the  model  “average  stay  time”  and  MJ^^  is  the  m^el  “travel  time”  for  each  type  unit. 
p,<4j  is  always  computed  from  the _NUFAM  III  input  values  and  M  ^j  rather  than  the  scenario 
average  stay  and  move  time  values  S  and  M. 


d.  If  we  enter  the  stay/move  process  at  some  random  time  T,  at  times  we  will  be  interested  in 
the  amount  of  time  remaining  after  T  until  the  unit  moves.  This  is,  of  course,  the  residual  life  of  the 
stay  time,  denoted  as  Y^.  The  distribution  of  is  given  in  Appendix  C. 

e.  When  the  POT  As  are  calculated  in  TADER,  it  relies  on  input  that  specifies  for  each  type  unit 
the  proportion  of  time  the  unit  is  stationary  and  moving  (broken  down  into  “in  the  open”  and 
“under  cover”).  Thus  the  POTA  is  already  weighted  for  the  relative  amount  of  time  any  given  type 
unit  spends  stationary  or  moving.  As  a  result,  we  do  not  check  in  NUFAM  III  to  see  whether  a  unit 
is  stationary  or  moving  at  the  time  that  an  acquisition  event  occurs.  The  proportion  of  time  spent 
stationary  or  moving  is  accounted  for  when  we  compute  the  probability  that  a  unit  is  available  for 
being  fired  upon  by  a  nuclear  weapon  (section  IV). 

f.  Example.  Suppose  that  a  unit  is  fairly  mobile  with  a  model  input  average  stay  time  S  of  15 

minutes  and  an  average  move  time  of  45  minutes.  Then  Pj^j,  =  =  0.25. 
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Section  II.  MODELING  THE  RETENTION  TIME 
2-3.  DEFINING  THE  RETENTION  TIME  T^ 

a.  The  retention  time  is  the  length  of  time  that  an  acquired  unit  will  remain  on  the  acquisition 
list.  A  fixed  tuget  remains  on  the  acquisition  list  indefinitely  (thus  the  methodology  discussed  in 
this  paper  does  not  apply  —  see  Youngren  (1988a]  ).  A  relocatable  target  unit,  on  the  other  hand, 
may  be  removed  from  the  acquisition  list  under  the  circumstances  discussed  below.  The  retention 
time  is  currently  represented  in  NUFAM  III  as  a  fixed  input  value  R.  R  (which  will  vary  by  tyjie 
unit  and  zone)  is  set  equal  to  E[T^]. 

b.  Although  all  relocatable  NUFAM  III  units  can  move  at  least  once  during  a  simulation,  some 
unit  types  are  engaged  while  moving  in  column.  We  refer  to  them  as  “units  engaged  while  moving”; 
this  designation  indicates  a  unit  in  march  column  detected  moving  as  part  of  a  longer  column  of 
similar  units.  We  refer  to  all  other  type  units  as  “engaged  while  stationary.” 

c.  Given  the  acquisition  systems  represented  in  the  TAS  III  study,  we  find  that  some  of  the 
short-range  sensors  which  can  reach  out  to  the  limits  of  Zone  3  (25  km)  can  provide  a  nearly 
continuous  tracking  capability  -  that  is,  they  can  continue  to  observe  a  unit  for  some  p)eriod  of  time 
after  acquisition.  The  tracking  may  be  continuous  or  nearly  continuous  (observed  frequently).  Units 
in  Zones  4  and  5,  on  the  other  hand,  cannot  be  tracked  after  acquisition.  The  distribution  used  to 
determine  the  retention  time  R  will  depend  upon  the  tracking  capability  and  whether  the  unit  is 
engaged  while  stationary  or  moving.  Three  cases  are  discussed  below. 

2-4.  TARGETS  ENGAGED  WHILE  STATIONARY  THAT  FACE  NO  TRACKING  CAPABILITY 

a.  This  case  models  the  situation  where  a  unit  is  acquired  while  stationary,  planned  for  fire  at 
that  location,  and  there  is  no  capability  to  track  the  location  of  a  unit  after  acquisition.  All  units  in 
Zones  4  and  5  except  those  which  are  identified  by  type  as  “engaged  while  moving”  fall  into  this 
category.  The  unit  may  be  acquired  while  stationary  or  moving;  however,  the  acquisition  is  useful 
for  fire  planning  only  if  it  was  acquired  while  it  was  stationary  or  when  it  was  moving  into  an 
identified  stationary  location  (e.g.,  assembly  area).  NUFAM  111  handles  this  by  checking  acquisitions 
of  units  of  this  type  to  see  if  it  was  stationary  at  the  time  of  the  last  observation,  denoted  by 
P/a#/  #/ar‘  units  that  were  stationary  at  the  time  of  the  last  observation  (in  this  case,  at  the 

time  of  acquisition)  are  planned  for  fire;  thus  the  retention  time  is  an  estimate  of  the  length  of  time 
that  the  target  unit  will  remain  at  the  location  at  which  it  was  acquired,  given  that  it  was  acquired 
while  stationary  or  when  moving  into  an  identified  stationary  location.  If  it  was  stationary  at  the 
time  of  acquisition,  this  will  be  the  residual  life  of  the  stay  time,  Y^.  If  the  unit  was  moving  into  the 
stationary  location  when  acquired,  it  will  have  the  entire  stay  period  S  ahead  of  it.  Thus 

E[Tfl  I  stay  ]  =  E[  Y^  |  stationary  at  time  of  acquisition  ] 

-I-  E[  S  1  moving  into  a  stationary  location  at  time  of  acquisition  while  moving  ] 

•  P[  moving  into  a  stationary  location  at  time  of  acquisition  |  moving  at  acquisition  ] 

=  E[Y5]  +  E[S]  •  6,  [6] 

where  6  represents  the  probability  that  a  unit  acquired  while  moving  is  moving  into  a  stationary 
location.  We  can  derive  6  from  the  scenario  stay/move  times  for  the  unit  by  considering  how  long  a 
period  is  possible  between  time  of  acquisition  of  a  moving  unit  and  the  time  it  stops  in  order  for  the 
stationary  location  to  be  known  at  the  time  of  acquisition.  We  denote  this  time  by  d\  thus  8  = 
P(  YJ^^  <  d  ],  where  Y^y^  is  the  residual  life  of  the  move  time.  Obviously,  d  will  be  small;  if  desired, 
d  can  be  approximated  as  zero. 


2-3 


CAA-TP-89-2 


b.  Since  d  is  assumed  to  be  small,  we  can  safely  assume  that  d  is  less  than  |  M  (  if  it  were  not,  a 
different  equation  for  P[  Ya/  <  d]  could  be  used;  the  equation  used  in  this  section  only  holds  for 
d  <  I  M  (see  Appendix  C)  ).  Thus,  P[  <  d  ]  =  =  (Appendix  C)  and 

E[Th  I  stay  ]  =  E[Ys]  +  E[S]  •  «  =  ^  S  +  S  •  1  =  S  •  [  ^  ^  .  [7] 

c.  Example.  Suppose  that  S  =  120  minutes  and  M  =  45  minutes.  If  d  =  5  minutes,  then  using 

the  formulas  in  Appendix  C,  P(  <  d  ]  =  5  min  /  45  min  =  0.11  and  E[  Y5  J  .=  (  14/27  )  • 

120  min  =  62.2  min.  Thus  E[  T^  |  stay  ]  =  120  min  •  (  0.52  +  0.11  )  =  75.6  min. 


2-5.  TARGETS  ENGAGED  WHILE  MOVING,  FACING  NO  TRACKING  CAPABILITY 

a.  This  case  represents  the  situation  where  a  unit  is  acquired  as  part  of  a  column  of  units 
moving  2dong  a  specified  route;  we  assume  no  capability  to  track  the  location  of  the  unit  after 
acquisition  (zones  4  and  5  using  TAS  III  data).  In  this  case,  the  nuclear  weapon  is  targeted  at  the 
location  of  the  unit  when  acquired  —  because  the  unit  is  one  of  many  moving  along  that  route,  it  is 
anticipated  that  a  similar  unit  will  be  within  the  radius  of  ejects  of  the  weapon  when  it  detonates. 
The  retention  time  is  an  estimate  of  the  length  of  time  that  other  units  will  be  moving  through  the 
location  at  which  the  unit  was  acquired.  An  estimate  of  this  time  in  hours  would  be: 

E[T^]  =  Expected  length  of  column  remaining  (km)  /  speed  of  column  (km/hr) 

=  (  Column  length  in  km  /  2  )  /  speed  of  column  (km/hr).  [8] 

b.  Example.  Suppose  that  a  column  moving  at  20  km/hr  is  2  km  long.  Then  E[  T^  ] 
=  [  (  2  km/2  )  /  20  km/hr  ]  =  0.05  hour  =  3  min. 

2-6.  TARGETS  FAQNG  A  NEARLY  CONTINUOUS  TRACKING  CAPABILITY 

a.  This  case  models  circumstances  where  there  is  a  capability  to  continuously  (or  intermittently 
with  a  short  duty  cycle)  monitor  or  track  the  location  of  an  acquired  unit.  All  units  placed  in 
NUFAM  III  zones  1,  2,  and  3  are  assumed  to  fall  into  this  category.  However,  this  tracking 
capability  cannot  last  forever  —  at  some  point,  the  unit  will  be  obscured,  a  gap  in  coverage  will 
occur,  etc.  In  this  case,  the  retention  time  is  an  estimate  of  how  long  the  unit  can  be  tracked  using 
the  same  sensors  that  initially  detected  the  unit.  Estimating  the  retention  time  under  these 
circumstances  is  discussed  below. 

b.  Approximating  a  continuous  or  nearly  continuous  observation  with  frequent  glimpses  is 
reasonable  if  the  time  between  glimpses  is  small.  All  units  in  NUFAM  III  are  acquired  using  the  2- 
hour  POTA  glimpse  cycle.  Thus  our  representation  of  tracking  must  use  the  same  2-hour  glimpse 
cycle.  Using  a  2-hour  glimpse  cycle  is  not  a  nearly  continuous  observation  model,  but  it  is  the  only 
model  that  can  be  based  on  current  POTA  data.  Ongoing  data  collection  efforts  at  CAA  (TAS  IV) 
will  permit  a  better  representation  of  continuous  observation  in  the  future. 

c.  A  unit  is  initially  detected  as  a  result  of  one  or  more  acquisition  cycles,  during  which  there  is 
a  probability  that  the  unit  will  be  detected  by  one  or  more  operational  sensors.  The  probability  of 
independently  detecting  any  given  unit  during  each  2-hour  cycle,  the  POTA,  is  implicitly  defined  as 
Prfdetect  during  this  cycle  (  no  detection  previous  cycle].  When  tracking  a  unit  after  it  has  been 
detected  using  the  same  sensors,  the  probability  that  at  least  one  of  the  sensors  will  continue  to 
detect  the  unit  is  Pr[detect  during  this  cycle  |  detection  previous  cycle],  which  would  be  a  type  of 
cued  detection.  We  denote  Pr[detect  during  this  cycle  ]  detection  previous  cycle]  as  Pq.  Clearly, 
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Pq  >  POT  A,  so  a  useful  representation  of  Pq  is  Pq  =  7  •  POTA,  where  1  <  7  <  1/POTA.  Thus, 
7  is  a  factor  that  indicates  how  much  better  the  sensors  are  at  detecting  a  unit  previously  detected 
and  identified  than  detecting  a  new,  previously  undetected  unit. 


d.  Just  as  the  use  of  the  POTA  in  NUFAM  III  leads  to  a  geometric  distribution  for  the  number 
of  cycles  before  a  target  is  acquired,  the  retention  time  should  have  the  same  distributional  form. 
The  probability  of  tracking  a  unit  for  k  acquisition  cycles  is  equal  to  (  7-  POTA  )  ^,  it  =  0,  1,  ...  , 
and  the  probability  that  a  unit  is  lost  (dropped  from  acquisition  list  )  on  the  (it  -f  1st)  cycle  is 
(  1  —  7  •  POTA  )  (  7*  POTA  )  ,  which  is  also  geometric.  Thus,  if  is  defined  as  a 

geometric  (  1  —  7  •  POTA  )  random  variable,  and  is  defined  as  a  uniform  (0,2)  random 
variable,  then  the  retention  time  T^  can  be  represented  as  2Xj^  +  U^.  Clearly 


EfTj  I  POTA  I  =  2  E(  X„  I  +  E[  1  =  , .^.^ppTA  )  +  ‘ 


1 

1  -  7-  POTA 


[9] 


e.  It  may  be  possible  to  estimate  Pq  (  i.e.,  7  •  POTA  )  by  running  the  TADER  model  using 
different  parametric  assumptions.  Alternatively,  we  handle  the  uncertainty  associated  with  the 
parameter  7  by  placing  a  prior  distribution  on  7  that  reflects  our  uncertainty.  Without  further 
information,  a  flat  prior  [  Uniform  (  1,  1/POTA  )  ]  is  reasonable.  Our  expected  value  for  Pq  using 
this  uniform  prior  halves  the  probability  that  a  subsequent  glimpse  will  fail  to  detect  the  target.  The 
expected  value  for  T^,  E[T^],  is  found  using  standard  techniques:  EfT^j]  =J-y  E[T^|7]  dG(7)  , 
where  G(y)  is  the  cumulative  distribution  function  (cdf)  of  7.  Unfortunately,  T^  does  not  have  a 
finite  mean  because  7  is  allowed  to  become  arbitrarily  close  to  1/POTA  (equivalent  to  Pq  =  1, 
which  yields  an  infinite  T^  ).  However,  NUFAM  III  is  intended  to  be  used  to  model  only  a  brief 
period  within  a  theater  battle  (not  to  exceed  12  hours).  With  the  24-hour  initial  acquisition  cycle 
used  in  NUFAM  III  runs  to  achieve  an  approximate  steady-state  acquisition  process,  the  simulation 
will  always  terminate  in  less  than  36  hours  of  simulated  time.  As  a  result,  any  time  T^  >  36  hours 
is  equivalent  to  an  arbitrarily  long  time  in  a  NUFAM  III  simulation  run.  Since  we  are  only 
interested  in  the  expectation  E[Tyj],  we  can  truncate  our  distribution  for  E[T^  I  7  ]  as  follows: 


E[Tfl  I  7,  POTA  1  _  1  _  y]  POTA  pAw  [  ^  A  ]  ’ 

=  36  for  7>p^[l-±].  [10] 

f.  We  calculate  that  [  1  -  J  =:  PC^  '  7  >  1,  7  can  be  <  only 

when  POTA  <  0.97;  thus,  E[Tfl  |  7  ]  =  36  when  POTA  >  0.97.  For  values  of  POTA  less  than 
0.97,  we  must  compute  the  expectation  of  T^: 

For  POTA  >  0.97,  E[Tfl]  =  36  hrs. 


For  POTA  <  0.97,  recalling  that  7  ~  Uniform  (  1,  1/POTA  ), 

0.97  1 

POTA  PoTa 

E[Tfll  =  1  -  7?  POTA  ■  {7=^% 


0.97 

P5TS 


=  1,= 


1 


pCta 


97  36 

P5Ta 


d7 


1 

P(5TA 


-  1 


7  =  1  1  -  POTA 


1 


POTA 


-7)-‘d7  ->-36 


1  -  0.97 

1  -  p6ta 
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After  some  calculation,  we  find  that 


EfT  1  —  POTA  —  0.97  I  r  1  —  0.97  ~[  ,  1-08  r_  pOTA  <  0  97’ 

E[Tj,]  =  36  for  POTA  >  0.97.  [11] 

We  use  this  value  as  the  retention  time,  R,  input  to  NUFAM  III  for  ail  units  in  Zones  1,  2,  and  3. 
g.  Example.  Suppose  that  a  unit  had  a  POTA  of  0.2.  Then 


0.2  - 

0.97  r  1 

-  0.97  1  . 

1.08  _  r  0.77  1 

(1  - 

0.2  )2  L  ] 

[  -  0.2  J 

1-0.2  “  L  6^  J 

Section  ffl.  COMPUTING  THE  PROBABILITY  THAT  AN  ACQUIRED  UNIT 
WILL  REMAIN  ON  THE  ACQUISITION  LIST 


2-7.  DEFINITION 


a.  The  probability  that  the  unit  on  the  acquisition  list  at  the  end  of  the  fire  planning  cycle  (time 
tp  as  shown  in  Figure  1-2)  will  remain  on  the  acquisition  list  until  some  elapsed  time  s  after  tp  is 
simply  the  survival  function  of  the  residual  life  of  the  retention  time  at  tp,  Y^,  evaluated  at  s. 
Recall  from  [2]  (paragraph  l-4b)  that 

P(  Yfl  >  s  ]  =  1  -  1  du  =  8.^  0  <  s  <  R;  0  for  s  >  R.  [12] 

Combining  this  probability  with  pae«  given  in  [2]  (paragraph  2- lb),  the  probability  that  a  unit  is  on 
the  acquisition  list  at  time  tp  and  remains  on  the  list  for  s  additional  hours  is 

P[  on  acquisition  list  at  tp  and  still  on  the  list  at  time  tp  s  ]  =  pacf  -  P[  Yy{  >  s  ] 


R 


POTA 


—  1  +  R 


R  —  a 


POTA 


-  1  R 


R  —  3 

~W~ 

0  <  a  <  R  ;  0  otherwise. 


[13] 


b.  In  NUFAM  III,  we  are  currently  using  a  1-hour  fire  execution  cycle,  with  targets  planned  for 
fire  if  they  are  acquired  at  any  time  up  to  the  beginning  of  the  fire  execution  cycle.  Thus,  the 
probability  that  we  acquire  a  unit  that  remains  available  for  engagement  through  the  end  of  a  1- 
hour  fire  execution  cycle  is 


^retain  7  —  acquired  and  is  retained  for  one  hour  ].  Setting  a  =  1  in  [13]  yields 


^retain  I 


R  -  1 


2 

P(^ 


-  1  +  R 


(  R>  1  ). 


c.  A  more  sophisticated  representation  recognizes  that  a  target  need  only  be  retained  until  it  is 
fired  upon.  The  time  of  detonation  is  denoted  as  Tj  (Figure  1-2).  Thus  the  remaining  (residual) 
retention  time,  Y^,  must  be  greater  than  Tj  — tp  in  order  for  the  target  to  be  fired  upon. 
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PI  >  T,  -  tp  I  T,  -  tp  ]  =  1  - 


—  tp< 


=  0 


for  Tj  -  tp  >  R. 


[14] 


Section  IV.  COMPUTING  THE  PROBABILITY  THAT  A  UNIT  IS  AVAILABLE  FOR  FIRE 
2-8.  DEFINITION 


a.  In  order  for  a  unit  to  be  available  for  engagement  (subject  to  fire)  by  nuclear  weapons,  it 
must  be  acquired  and  retained  as  a  target  at  least  until  the  scheduled  firing  time.  The  target  is 
ev2duated  at  some  known  planning  time  tp;  during  the  fire  planning  cycle,  the  detonation  will  be 
scheduled  to  occur  at  some  time  T^  within  the  fire  execution  cycle.  Thus  a  target  that  is  on  the 
acquisition  list  at  time  tp  must  be  retained  at  least  until  T^;  this  occurs  when  the  residual  time  on 
the  acquisition  list,  Y^,  is  greater  than  or  equal  to  (  T,^— tp  ).  We  define  the  probability  P[  unit  is 
acquired  and  retained  until  detonation]  as  Thus 

Yfl  >  T<j— tp  I  acquired  at  tp  ]  •  P[  acquired  at  tp  ]. 


b.  In  the  absence  of  further  information  about  the  firing  times,  it  is  reasonable  to  assume  that 
any  detonation  time  within  the  fire  execution  cycle  is  equally  likely.  If  the  fire  execution  cycle  is  1 
hour  long,  and  the  cycle  begins  immediately  after  the  planning  cycle  ends,  then  (T,^— tp)  ~  U[  0,1  ]. 
If  a  different  distribution  for  (T<j— tp)  can  be  inferred  or  if  the  time  at  which  the  fire  planning  cycle 
ends  is  not  the  time  at  which  the  fire  execution  cycle  begins,  the  probability  Y/j  >  (Tj  —  tp)  is 
computed  in  the  same  manner  as  shown,  but  a  different  set  of  limits  on  the  integral  and  a  different 
formula  for  P[Tj— tp  =  t]  will  be  required.  We  calculate  using  the  residual  life  of  the 

retention  time  : 


P[  Y^  >  T,j  —  tp  ]  —  P[  Yyj  >  Tj  —  tp  ]  P[  T^  —  tp  —  t  ]  dt  — 

where  D  =  min  {  R,  1  }  . 

Thus 

P(Y„>T,-t,I  =  |J’  =d-^  =  D(I.  5%). 

Therefore, 


P  avaiV 


~T 
pDta 


R 


-  1  +  R 


[15] 


This  yields: 

^  1=  Pavo.-/ 

R  ^  P  avail 


R^ 


pSta 


-  1  -I-  R 


2R-  1 


P(5TA 


-  1  +  R 


[16] 


c.  In  some  cases,  all  weapons  of  a  particular  type  will  detonate  before  the  end  of  the  fire 
execution  cycle.  An  alternative  distribution  for  (T^— tp)  that  requires  no  additional  calculation  is 
(Tj  — tp)  ~  uniform  [  0  ,  L  ],  where  L  is  the  time  to  the  last  detonation  after  the  beginning  of  the 
fire  execution  cycle  (  L  <  1  hour).  Thus  D  =  min  {  R,  L  }  and  equation  [15]  can  be  used  to 
calculate  This  approach  may  also  be  used  to  generalize  the  results  for  fire  execution  cycles  of 

length  L,  where  L  ^  1  hour. 
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d.  Example.  Suppose  that  the  POTA  is  0.2  and  the  retention  time  R  is  5.3  hours.  Using  a  1- 
hour  fire  execution  cycle,  D  =  min{  1  hour,  5.3  hours  }  =  1  hour  and 


P  avail 


2R-  1 


2(5.3)  -  1 


2 

P(5TA 


—  1  +  R 


T2 

'Lo:5  " 


1  +  5.3 


Section  V.  COMPUTING  THE  PROBABILITY  THAT  A  UNIT  HAS  MOVED 
AT  THE  TIME  OF  DETONATION 

2-9.  DEFINITION.  Because  we  cannot  track  acquired  units  after  a  certain  point  (that  point 
depending  upon  our  capability  for  near  continuous  tracking),  it  is  possible  for  a  weapon  to  miss  a 
target.  We  plan  the  weapon  to  detonate  at  the  Desired  Ground  Zero  (DGZ)  at  which  the  weapon 
effects  against  the  target  are  maximized  (subject  to  some  physical  constraints).  If  a  target  unit  has 
moved  since  the  time  the  last  DGZ  adjustment  was  completed,  the  weapon  will  miss  the  target.  Unit 
movement  is  represented  through  the  NUFAM  III  inputs  “average  stay  time”  and  the  “travel 
time”  M discussed  previously.  The  values  for  S and  M  are  determined  differently  for  different 
tracking  capabilities. 

2-10.  TARGETS  ENGAGED  WHILE  STATIONARY,  FACING  NO  TRACKING  CAPABILITY. 


a.  When  we  have  no  capability  to  track  units,  we  can  only  successfully  engage  them  if  they  are 
stationary  at  the  time  of  acquisition  (which  is  assumed  to  hold  for  units  in  this  category  by 
definition)  and  they  remain  in  place  at  the  point  at  which  they  were  acquired.  The  average  sUy  time 
input  to  the  model  equals  the  scenario  average  time  the  unit  will  remain  stationary  (S); 
equals  the  scenario  average  time  that  a  unit  will  be  moving  from  stationary  location  to  stationary 
location  (M). 

b.  NUFAM  III  will  check  the  probability  that  the  unit  was  stationary  at  the  time  it  was  last 
observed  before  launching  the  weapon.  We  denote  this  probability  as  Pi^,f  it  is  stored  in  the 
TUNIT  file  as  PROB.STAY.LOBS.  The  last  observation  for  units  in  this  category  was  at  the  time 
of  acquisition,  so  referring  back  to  paragraph  2-4a, 


PROB.STAY.LOBS  =  =  P[  stationary  at  time  of  acquisition  ] 

-i-  I  P[  moving  into  a  stationary  location  at  time  of  acquisition  |  moving  at  acquisition  ] 

•  P[  moving  at  acquisition  ]  | 

[17] 


_  S  _ 
S  -H  M 


-h  6 


M 

S  M’ 


where  6  =  P[  Yj^  <(/]  =  i,  i  <  |m  (Appendix  C). 

c.  The  probability  that  the  unit  remains  stationary  at  its  identified  location  until  the  time  of 
detonation  is  equal  to  the  probability  that  the  time  that  the  unit  remains  at  the  location  at  which  it 
was  acquired  is  greater  ^han  the  time  from  acquisition  to  detonation.  If  we  denote  the  time  of 
acquisition  as  Tacf  and  the  time  of  detonation  as  T^,  then  the  probability  that  the  unit  remains  at 
its  identified  location  until  the  time  of  detonation  is  P(  Y5  >  Tj  —  Tacf  and  stationary  at  Tac^  ]. 
The  time  Tac^  is  not  known,  but  the  difference  (  tp  —  Taej  )  is  distributed  as  the  age  of  the 
distribution  of  the  time  of  acquisition,  denoted  as  Y^j^,  independent  of  the  distribution  of  T^  — tp. 
The  variable  Tj  —  tp  can  have  any  nonnegative  distribution;  this  paper  uses  either  a  constant  (rj)  or 
uniform  (  0,  L  )  distribution,  where  L  is  the  length  of  the  fire  execution  cycle  and  tp  is  the  time  that 
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the  fire  planning  cycle  ends  (and  the  fire  execution  cycle  begins).  The  formula  for  determining 
P[  Y5  >  Tj  —  Tacj  and  stationary  at  Tac»  li  taking  into  account  the  distributions  of  (  —  tp  ) 

and  (  tp  —  Tacj  ),  may  be  found  in  Appendix  C. 


d.  To  simplify  matters,  the  current  studies  using  NUFAM  III  determine  the  probability  that  the 
unit  is  stationary  at  that  identified  location  at  the  time  of  detonation  by  drawing  against  p^^^y  (  [5], 
paragraph  2-2c  ).  This  is  equivalent  to  ignoring  the  move/stay  state  at  the  time  of  acquisition  Tacf 
or  any  time  between  Tccf  and  T^.  Given  the  definition  of  these  units  (units  acquired  and  engaged 
while  stationary),  this  is  not  unreasonable.  Thus 

.<«,  =  1*  P.ia,  =  Sat  =  S,  and  =  M.  [18] 

Future  studies  may  set  ,tay  ~  g — +M —  ~ 

P[  Y5  >  Tj  —  Tflcj  and  stationary  at  T»cj  ]. 

2-11.  TARGETS  ENGAGED  WHILE  MOVING,  FACING  NO  TRACKING  CAPABILITY 


a.  When  we  have  no  capability  to  track  units  detected  while  moving,  we  cannot  successfully 
engage  them— however,  if  the  acquired  unit  is  part  of  a  column  moving  along  a  defined  path,  we  may 
be  able  to  successfully  engage  a  different  unit  that  subsequently  moves  through  the  point  at  which 
the  original  unit  was  acquired.  The  average  stay  time  S input  to  the  model  represents  the  average 
time  that  a  unit  is  passing  through  the  acquisition  point;  represents  the  average  time  that  no 
unit  is  passing  through  the  acquisition  point.  If  the  column  is  composed  of  several  units  of  the  same 
type  (the  basic  assumption  defining  a  type  unit  “engaged  while  moving”),  then  we  may  compute 
Sa/  and  Si  A/  as  follows: 


Sa/ 


length  of  each  unit  in  column  (km) 
speed  of  column  (km/hr) 


length  of  gaps  between  units  in  the  column  (km) 
speed  of  column  (km/hr) 


[19] 


b.  The  probability  that  the  unit  was  stationary  at  the  time  it  was  last  observed  immediately 

before  firing  the  weapon,  ,  is  not  applicable  for  units  in  this  category.  As  a  result,  we  set 

PROB.STAY.LOBS  =  p,„,  =  1. 

c.  Example.  Suppose  that  the  length  of  each  unit  is  250  meters  =  0.25  km,  and  the  column  is 
traveling  at  20  km/hr.  Suppose  that  the  length  of  the  gaps  between  each  unit  is  50  meters  =  0.05 
km.  Then  Sat  =  0.25  /  20  =  0.0125  hr.  and  =  0.05  /  20  =  0.0025  hr. 

2-12.  TARGETS  FACING  A  NEARLY  CONTINUOUS  TRACKING  CAPABILITY 
(SCHEDULED  FIRE) 

a.  We  have  assumed  that  the  target  units  of  this  type  may  be  tracked  until  shortly  before  the 
weapon  detonates  at  the  time  T^;  Tj  having  been  scheduled  (set  during  the  time  of  fire  planning). 
The  tracking  is  discontinued  shortly  before  firing  due  to  the  necessity  of  communicating  the  fire 
order  to  the  firing  unit  and  adjusting  the  delivery  system  for  fire.  Thus  the  DGZ  of  the  weapon  has 
been  shifted  along  with  the  unit  until  the  order  to  fire  is  given.  From  another  perspective,  the  unit 
(whether  stationary  or  moving)  has  remained  stationary  with  respect  to  the  DGZ  during  this  entire 
tracking  period.  As  a  result,  the  ratio  Sa^  /  [Saj  +  Ma^]  denotes  the  probability  that  the  unit 
remains  stationary  with  respect  to  the  last  computed  DGZ.  If  it  does,  then  it  will  be  subject  to 
nuclear  effects.  If  not,  the  target  unit  must  have  shifted  unexpectedly  during  the  interval  between 
final  fire  adjustment  and  firing,  causing  the  weapon  to  miss. 
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b.  Let  r}  denote  the  amount  of  time  between  the  final  observation  (time  of  final  DGZ 
adjustment)  and  the  detonation.  This  delay  accounts  for  the  time  necessary  to  transmit  the  order  to 
fire  to  the  firing  unit,  adjust  the  delivery  system,  etc.  If  the  unit  is  stationary  at  the  time  that  the 
last  observation  is  made,  then  the  DGZ  will  be  the  final  location.  The  weapon  will  miss  if  the  unit 
moves  during  the  period  rj  (Note:  If  the  unit  moves  at  a  time  slightly  less  than  rj,  it  may  still  be 
within  the  effects  radii  for  the  burst,  but  this  adjustment  is  dependent  on  weapon  as  well  as  target 
characteristics;  its  contribution  to  the  probability  is  normally  small  for  small  yield  weafwns  and  will 
be  ignored).  Thus  the  probability  that  the  unit  has  not  moved  from  the  location  at  which  it  was  last 
observed,  denoted  by  pno  movci  is  P[  Y5  >  tf  and  the  unit  stationary  at  time  of  last  observation  ], 
which  is: 


Pno  move  =  P[  Yj  >  i;  |  stationary  at  time  of  last  observation  ] 


c.  Let  S  and  M  denote  the  scenario  average  stay  and  move  times  for  a  particular  type  of  unit 

2  c 


under  consideration.  For  ij  <  |  S  (Appendix  C), 

P[  Y^  >  »j  I  stationary  at  time  of  last  observation  ]  =  !—».  Thus 

3 


Pno  moot 


= 

S 


_  S  _  _  _S-T?_ 
S  +  M  S  +  M 


[20] 


d.  NUFAM  III  will  check  the  probability  that  the  unit  was  stationary  at  a  time  t}  before 
detonation  during  the  fire  execution  cycle  immediately  before  firing  the  weapon.  We  denote  this 
probability  as  computed  as  follows: 


S 


PROB.STAY.LOBS  =  P[  stationary  at  last  observation  ]  =  pj^^,  = 

5 


[21] 


e.  When  NUFAM  III  computes  the  probability  that  a  detonated  round  missed  the  target  unit,  it 
checks  to  see  if  the  unit  has  moved.  Thus  to  force  the  current  NUFAM  III  code  to  record  a  miss  with 
the  appropriate  probability,  the  values  and  Mjj^  entered  into  the  data  base  must  solve 


_ 


=  P[  Y5  >  7  I  stationary  at  time  of  last  observation  ] 


Pno  move  _  S  —  Tf 

P/oit  tiay  S 


We  suggest  setting  =  rj  and  =  S  —  7. 

f.  Example.  Suppose  that  S  =  15  min,  M=  45  min,  and  7  =  10  min,  which  is  <  |  S, 
thus  P[  Y^  >  7  ]  “  1  “*  ^  and  pno  move  *  —  0.083.  Sj^  —  S  "“J/  —  5  min  and 

MJ^£  =  7  =  10  min. 


2-13.  TARGETS  FACING  A  NEARLY  CONTINUOUS  TRACKING  CAPABILITY 
(WITHOUT  SCHEDULED  FIRE) 

a.  This  case  represents  the  same  situation  as  paragraph  2-12  where  at  least  one  sensor  can 
observe  or  “track”  the  target  continuously  (or  nearly  continuously)  after  acquisition.  However,  it  is 
reasonable  to  suppose  that  if  we  can  track  the  target  unit  continuously,  we  would  give  the  order  to 
fire  when  the  unit  stops  for  the  first  lime  within  the  fire  execution  period.  Until  now,  we  have  given 
the  order  to  fire  (determined  the  time  Tj)  at  time  tp.  In  this  case,  we  define  a  (random)  time  to 
give  the  fire  order,  T y,  during  the  fire  execution  period  which  depends  upon  the  time  that  the  unit 
stops.  It  is  reasonable  to  assume  that  some  time  will  need  to  be  spent  getting  the  order  to  the 
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delivery  unit,  adjusting  the  aimpoint  to  the  location  at  which  the  unit  stopped,  confirming  that  the 
location  meets  the  criteria  for  employment  of  the  weapon,  etc.  As  in  the  previous  section,  this  time 
(T^— Ty)  is  denoted  as  7,  and  for  simplicity,  Tf  is  assumed  constant.  If  the  length  of  the  fire 
execution  cycle  is  L,  then  clearly  0  <  »j  <  L. 

b.  If  the  unit  is  stationary  at  time  tp,  T y=tp;  otherwise,  T y  equals  the  first  time  the  unit  stops 
between  tp  and  tp+(L— t;)  (Figure  2-2).  If  the  move  time  distribution  is  not  bounded  from  above  at 
some  value  less  than  the  fire  execution  period  duration,  there  is  a  positive  probability  that  the  unit 
will  be  moving  during  the  entire  fire  execution  period.  Therefore,  the  stay  and  move  times  Sj^  and 
M  ^  must  be  actjusted  to  account  for  this.  Because  the  stay/move  process  is  assumed  independent  of 
the  acquisition  process,  Ty  is  random  with  respect  to  the  acquisition  process.  Therefore,  the 
probability  that  the  unit  is  on  the  acquisition  list  at  the  time  Ty  is  simply  Paef 


Acquired 

Not  Acquired  T  4 

_Li_ 

Stay  S 

s 

Move  M 

M 

-  L  - 

1  Fire  Planning  Cycle  ] 

Fire  Execution  Cycle  1 

1 

t«f-q  — ' 

Time 

Tpef  tp 

T;  T, 

Time 


Figure  2-2.  Engaging  a  Target  Unit  Faced  with  Nearly  Continuous  Tracking,  without  Scheduled  Fire 

c.  The  probability  that  the  unit  is  available  for  fire,  given  that  it  was  on  the  acquisition  list  at 

Ty  is  still  the  probability  that  the  acquisition  is  stationary  at  Ty  and  retained  between  Ty  and  T^. 
However,  if  we  issue  the  fire  order  at  time  Ty,  we  will  retain  the  unit  between  Ty  and  T^  by 
definition.  Thus  the  conditional  probability  P[  unit  retained  until  T^  |  unit  on  the  acquisition  list  at 
Ty  ]  =  1.  Thus  =  1  •  p.c,.  [22] 

d.  Let  T]  denote  the  amount  of  time  between  the  final  observation  (time  of  final  DGZ 
adjustment)  and  the  detonation.  Since  the  unit  is  stationary  at  the  time  that  the  last  observation  is 
made  by  definition,  then  the  DGZ  will  be  the  final  location.  The  weapon  will  miss  if  the  unit  moves 
during  the  period  q.  If  Ty  =  tp,  the  unit  was  somewhere  in  its  stay  period  at  time  tp  and  the 
probability  that  it  will  remain  in  place  until  Tj  is  P[  Y5  >  q  ],  where  Yj  is  the  residual  life  of  the 
stay  time.  If  tp  <  Ty  <  tp+(L  — q),  then  the  unit  stopped  during  the  fire  execution  period  and  the 
time  Ty  is  equal  to  the  time  that  the  unit's  stay  period  begins.  Thus 
P[  stationary  atT^  |  tp  <  Ty<  tp-t-(L  — q)]  =  P[S  >  qj. 
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e.  The  probability  that  the  unit  was  stationary  at  time  tp  or  stopped  during  the  fire  execution 
cycle,  and  it  has  not  moved  from  the  location  at  which  it  was  last  observed  (pno  move)  is: 


Pno  mo»e=  |  P[  Y5  >  »j  |  stationary  at  tp  ] 


S  +  M 


+  P[  S  >  T/  1  stopped  at  T moving  at  tp  ]  •  P[  Ym  <  L  — fj  |  moving  at  tp  ]  •  -  ^  | 

S  -f-  M  J 


=  P[  Y5  >  »?  ] 


S  +  -M 


+  P(S  >  f?I  •  P(  Ym  <  L-fj] 


M 


S  +  M’ 


[23] 


P/a*<  probability  that  the  unit  was  stationary  r)  time  units  before  detonation  during 

the  fire  execution  cycle  immediately  before  firing  the  weapon,  is: 

PROB.STAY.LOBS  =  +  P[  Y,„  <  L-,,  ]  •  [24] 

g.  When  NUFAM  III  computes  the  probability  that  a  detonated  round  missed  the  target  unit,  it 
checks  to  see  if  the  unit  has  moved.  Thus  to  force  the  current  NUFAM  III  code  to  record  a  miss  with 
the  appropriate  probability,  the  values  S ^  and  M entered  tnio  the  data  base  must  solve 


_  Pno  move 

ttaf 

_  P(  Y5  >  r,  1  .  +  P[  S  >  r,  ]  .  P[  Yo,  <  L-i,  ]  . 


M 


S  +  M  '  S  +  M 

We  suggest  setting  S^f  =  ?[  >  r,  ]  •  ^  +  P[  S  >  »;  ]  •  P[  Ym  <  L-rj  ]  •  - - = 

b  +  M  S  -f  M 

Pno  mo»e  and  =  P[  Y5  <  rj  ]  •  J  ^  -  +  P[  S  <  rj  ]  •  P[  Ym  <  L  — 17  ]  •  55 

o  +  m  o  +  .M 


P/OJI  fiat  ~  p»<>  ”“»»«• 


h.  Example.  Suppose  7  =  10  min,  L  =  1  hour.  S  =  15  min,  and  M  =  45  min.  Then  (  L  — i;  ) 
=  50  min.  From  Appendix  C, 


P[  Ym  <  L-r;  ]  =  1  -  P[  Ym  >  L-r;  ] 


=  j^2(L-T7)  -  — ~  3*  30  min  -  |  M  ^  (  L-17)  <  |  M  =  60  min, 

=  0.962.  Also  -  ■-  =  0.25  and  -  ^  =  1  -  0.25  =  0.75 

S  +  M  S  4-  M 


Thus  p,„,  ^ 


f.  +  P[  Ym  <  L-r?  ]  •  y  ^  =  0.25  +  0.962  •  0.75  =  0.972. 

S  +  M  S  -t-  M 

Pno  move=  P(  Y5>r?  ]  •  »— §-=  +  P[  S  >  r?  ]  •  P[  Ym  <  L  — T)  )  •  >  -a 

=  ^  •  0.25  +  1  •  0.962  ■  0.75  =  0.805. 

S^lf  =  Pno  move  =  0.805,  M  —  Pno  move  =  0.972  -  0.805  =  0.167. 
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Section  VI.  COMPUTING  THE  PROBABILITY  THAT  A  UNIT 
IS  FIRED  AT  BY  A  NUCLEAR  WEAPON 


2-14.  DEFINITION 

a.  The  probability  that  a  unit  can  be  fired  at  by  a  nuclear  weapon,  denoted  as  is  the  joint 

probability  that  a  unit  is  on  the  acquisition  list,  is  not  dropped  from  the  acquisition  list  prior  to 
detonation,  and  is  stationary  at  the  time  of  last  observation.  This  probability  calculation  implicitly 
assumes  that  a  weapon  of  the  appropriate  yield  is  available  and  will  be  fired  at  the  target  if  it  meets 
these  criteria.  Recalling  that  =  P[  unit  is  acquired  and  retained  until  detonation]  , 


P fired  ^ avail  P/aj<  tiay 

P/irei~P“«»‘  Iq  >  tj— tp  I  Tj  — tp  —  t  ]  P[  Tj  — tp  =  t  ]  dt  ■  Pinai  ,iay 

Units  in  Zones  4  and  5,  regardless  of  whether  they  are  “engaged  while  stationary”  or  “engaged  while 
moving”,  are  assumed  to  be  stationary  or  moving  by  definition  at  the  time  of  acquisition,  and  they 
are  not  tracked  during  the  interval  between  acquisition  and  firing.  Thus  the  probability  that  they 
where  stationary  with  regard  to  the  location  at  which  they  were  ^u:quired,  can  be  ignored, 

so  it  is  set  equal  to  1.  Units  in  Zones  1,  2,  and  3  are  tracked  during  the  interval  between  acquisition 
and  firing,  thus  the  probability  p,^^^  is  equal  to  p^j^y  (scheduled  fire)  or  is  equal  to  + 

P[  Ym  <  L  — ]  •  pmove  when  the  fire  order  is  not  planned. 


b.  For  units  in  Zones  1,  2,  and  3,  scheduled  fire,  from  [15]  and  [21]  we  get 


^ fired 


R 


pSta 


-  1  -1-  R 


D  1  1 


2R'  S  +  M’ 


for  units  in  Zones  1,  2,  and  3,  without  scheduled  fire,  from  [15]  and  [24]  we  get 


P/iVei  —  5 

p?5ta 


- -  .  D  [  1  -  A  ]  .  r  -h  P[  Y„,  <  L-,,  ]  .  , 

_1  +  R  S+M_ 


and  for  all  units  in  Zones  4  and  5,  from  [15]  and  [18]  we  get 
-  D  [  1  -  ;^  ]  •  1  , 


^ fired 


pCta 


-  1  R 


[25] 


[26] 

[27] 


where  D  =  min{  R,  L  },  L  is  the  length  of  the  fire  e.xecution  time  (hours),  R  is  the  retention  time  in 
hours,  T)  is  the  length  of  time  between  the  time  the  fire  order  is  given  (time  the  unit  was  last 
observed  when  tracked)  and  the  time  of  detonation,  and  S  and  M  are  the  scenario  average  stay  and 
move  times  respectively. 


c.  Example.  Suppose  ?;  =  10  min,  L  =  1  hour.  S  =  15  min.  and  M  =  45  min  as  given  in  the 
previous  example.  Also  suppose  that  the  POTA  is  0.2  and  R  =  5.3  hours.  Then  for  units  in  Zones  1. 
2.  and  3  with  scheduled  fire. 


^ fired 


5.3 


-2-  -  1 

0.2 


+  5.3 


1  f  1  - 


15 


2- 5.3 


154-  45 


=  0.371  •  0.906  •  0.25  =  0.084. 
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For  units  in  Zones  1,  2,  and  3  without  scheduled  fire, 

=  X-l+5.3  ■  ‘  •  isfis]  =  ■  »•“«  ■ 

For  units  in  Zones  4  and  5  using  the  same  S  and  M  values  as  used  in  Zones  1,  2,  and  3, 

t  ■  '  t  ‘  ■  jij I  ■  ‘  = 

jj:5  -  1  +  5.3 


Section  VII.  COMPUTING  THE  PROBABILITY  THAT  A  UNIT 
IS  HIT  BY  A  NUCLEAR  WEAPON 

2-15.  DEFINITION 

a.  The  probability  that  a  unit  is  hit,  denoted  as  ,  is  the  joint  probability  that  a  unit  is  on 
the  acquisition  list,  is  not  dropped  from  the  acquisition  list  prior  to  detonation,  is  stationary  at  the 
time  of  last  observation,  and  does  not  move  from  the  last  computed  OGZ.  This  calculation 
implicitly  assumes  that  a  weapon  is  available  for  firing  at  the  unit,  the  weapon  is  fired,  and  both  the 
weapon  system  and  warhead  function  correctly.  Thus 


PAi<  —  Pa»ai7  ‘  P«» 

P\i«  =  P*'?  ■  Iq  I  tp  =  t  ]  P[  Tj  — tp  =  t  )  dt  •  p,^^, 

D  =  D  • 

PhU  PfireJ  c  4-  M  ’ 


Pno  move 


[28] 


where  Py,„^  is  computed  as  given  in  the  previous  section  and  SJ^f  and  are  the  values  entered 
into  the  NUFAM  III  data  base  for  the  unit  stay  time  and  move  time  respectively. 


b.  Example.  Suppose  the  probabilities  that  a  unit  will  be  fired  on  are  as  given  in  the  previous 
example,  with  and  as  calculated  in  the  examples  at  the  end  of  each  paragraph  in  Section  V. 
For  the  purpose  of  comparison,  the  same  M  and  S  values  are  used  for  all  of  the  e.xamples  below.  In 
actual  scenarios,  the  M  and  S  values  would  vary  considerably  by  zone,  type  unit,  and  unit  activity. 


(1)  For  units  in  Zones  I,  2,  and  3  with  scheduled  fire,  p^„,f  =  0.084,  =  »7  =  10  min, 

S^=S  — >7  =  5  min,  and  p^,^  =  0.084  •  (  5/15  )  =  0.028. 

(2)  For  units  in  Zones  1,  2,  and  3  without  scheduled  fire,  =  0.319,  Sj^f  =  Pno  move  = 

0.805,  -  Pno  move  =  0.167,  and  p^-,  =  0.319  •  (  0.805/0.972)  =  0.264. 

(3)  For  units  in  Zones  4  and  5  engaged  while  stationary,  Pfi^ed  ~  0.336,  M  =  M  =  45 
min,  S  =  S  =  15  min,  and  pjj,y  =  0.336  •  (  15/60  )  =  0.084. 

(4)  For  units  in  Zones  4  and  5  engaged  while  moving,  ~  0.336.  S^f  =  0.0125  hr.  M 

=  0.0025  hr.  and  p^,,  =  0.336-  (  0.0125/0.0150  )  =  0.280. 
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CHAPTER  3 
SUMMARY 


SUMMARY. 

a.  The  principal  findings  of  the  work  reported  in  this  paper  are: 

(1)  The  target  acquisition  and  unit  movement  processes  may  be  represented  using 
alternating  renewal  processes. 

(2)  The  current  data  for  target  acquisition  used  within  NUFAM  III,  the  probability  of  target 
acquisition  (POTA)  values  for  a  2-hour  glimpse  period  derived  from  the  Target  Acquisition  Study, 
Version  III  (TAS  III),  require  the  number  of  periods  until  acquisition  to  be  distributed  as  geometric. 

(3)  The  process  of  tracking  nuclear  targets  after  acquisition  may  be  represented  based  on 
POTA  values  if  certain  assumptions  are  accepted. 

b.  The  efficiency  of  NUFAM  III  can  be  increased  by  eliminating  the  tedious,  time  consuming 
simulation  of  target  acquisition  cycles  by  drawing  against  probabilities  of  acquisition,  engagement, 
hit,  etc.  This  representation  allows  us  to  easily  compute  offline  the  distribution  of  the  number  of 
units  of  a  given  type  that  will  be  available  for  nuclear  engagement  for  any  pair  of  POTA  and 
retention  time.  We  can  use  this  to  determine  the  effect  of  variations  in  the  POTA  and  the  retention 
time  on  the  acquisition,  engagement,  and  defeat  of  various  units  without  running  NUFAM  III. 

c.  The  accuracy  of  the  NUFAM  III  output  can  be  enhanced  when  the  input  values  S,  M  and  R 
(the  stay  time,  move  time  and  retention  time  respectively)  are  computed  for  each  type  unit  using 
the  equations  documented  in  this  paper.  This  permits  the  NUFAM  III  model  to  distinguish  between 
different  types  of  target  units  and  different  sensor  capabilities,  yielding  a  correct  representation  of 
the  probability  that  a  unit  can  be  engaged  with  a  nuclear  weapon. 

d.  This  technique  may  also  be  applied  to  support  the  nuclear  effects  module  for  the 
IWFORCEM  model,  NEMESIS,  until  more  accurate  data  and  methodology  are  available.  We  can 
compute  the  acquired  target  sets  that  are  eligible  for  engagement  by  computing  Pj„j,/  as  stated 
above.  The  probability  that  acquired  units  are  hit  may  be  calculated  using 
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APPENDIX  B 

GUIDANCE  TO  NUFAM  III  USERS 

Section  I.  PARAMETERS 


B-1.  INTRODUCTION 

In  this  appendix,  we  discuss  how  to  determine  the  values  for  the  average  stay  time,  the  travel 
time,  and  the  retention  time  for  each  type  unit  input  to  NUFAM  III. 

B-2.  PHYSICAL  PARAMETERS  PROVIDED  BY  DOCTRINE  OR  THE  SCENARIO 

S  =  the  scenario  average  time  the  type  unit  will  remain  stationary. 

M  =  the  scenario  average  time  that  a  unit  will  be  moving  from  stationary  point  to  the  next 

stationary  point. 

8  =  the  probability  that  a  unit  detected  while  moving  is  moving  into  a  stationary  location. 

i  =  how  soon  a  unit  acquired  must  arrive  at  its  stationary  location  in  order  for  the 
stationary  location  to  be  known,  if  the  unit  is  detected  while  moving  (i.e.,  the 
maximum  interval  of  time  between  detection  and  arrival  for  the  arrival  location  to  be 
identified). 

r\  =  the  amount  of  time  between  the  final  observation  (  time  of  final  fire  adjustment  )  and 
the  detonation  when  tracking  a  target  unit.  This  is  primarily  a  function  of  the 
command  and  control  time  associated  with  firing  a  nuclear  round  from  the  designated 
firing  unit. 

7  =  a  multiplicative  factor  that  indicates  how  much  better  the  sensors  are,  relative  to  the 

POTA,  at  continuing  to  track  a  unit  previously  detected  and  identified  than  detecting 
a  new,  previously  undetected  unit. 

L  =  the  length  of  the  fire  execution  cycle. 

B-3.  NUFAM  III  INPUT  VALUES 

Sjj^  =  AVERAGE.STAY.TIME  (by  type  unit  and  zone) 

MJ^^=  TRAVEL.TIME  (by  type  unit) 

R  =  RETENTION.TIME  (by  type  unit) 

NOTE:  Parameters  in  paragraph  B-2  are  stated  in  terms  of  hours,  which  lead  to  the  calculation  of 
S]^f,  and  R  in  hours.  Before  entry  into  the  NUFAM  Ilf  input  data  base,  Sj^f,  and  R  must 
be  converted  to  minutes  (by  multiplying  by  60).  Also  note  that  NUFAM  III  will  compute  S  and  M 
for  each  unit  in  minutes,  using  S^and  Mj^.  The  actual  values  of  and  Mj^  used  are  not  critical; 
what  matters  is  tlw  ratio  Sjj£  /  (  -b  Mjj^  ).  Therefore,  in  order  to  avoid  roundoff_probletm,  the 
computed  values  and  My^  should  be  multiplied  by  some  constant  such  that  +  Mj^  is 
approximately  equal  to  120  min. 
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Section  H.  MAJOR  ASSUMPTIONS 


B-4.  ASSUMPTIONS  MADE  WITHIN  THE  NUFAM  HI  MODEL 

a.  POTA  values  from  TAS  III  represent  the  independent  probability  of  acquisition  during 
successive,  2-hour  glimpse  periods. 

b.  Units  are  acquired  and  move  independently  of  each  other. 

c.  The  end  of  the  fire  planning  time,  tp,  is  set  independently  of  the  acquisition  and  move/stay 
processes.  The  fire  execution  cycle  begins  at  tp  and  has  a  length  of  L  (  L  =  1  hour  in  current 
studies  using  NUFAM  III). 

d.  The  number  of  2-hour  periods  required  to  detect  a  unit,  is  distributed  as  Geometric(POTA). 

e.  For  a  detected  unit,  the  specific  time  of  acquisition  in  the  (2-hour)  period  detected,  is 
distributed  as  Uniform  (0,2). 

f.  Retention  time,  R,  is  set  to  a  constant  value  equal  to  E(T|^)  =  average  retention  time. 

g.  All  sensor  suites  observing  target  units  in  zones  1-3  have  nearly  continuous  tracking 
capability. 

h.  All  sensor  suites  observing  targets  in  zones  4  and  5  are  incapable  of  tracking  after 
acquisition. 

i.  Within  the  model,  the  duration  of  stay  time  for  a  unit  is  distributed  as 

Uniform(  ),  where  Sjy/  =  average  stay  time  input  to  NUFAM.  The  scenario  stay  time  S 

is  assumed  to  be  distributed  as  Uniform(  |§,  )  where  S  =  average  stay  time  for  the  type  unit  in 

the  scenario. 

j.  Within  the  model,  the  duration  of  move  time  for  a  unit  is  distributed  as 

Uniform(  ),  where  ”  average  move  time  input  to  NUFAM.  The  scenario  move 

time  M  is  assumed  to  be  distributed  as  Uniform(  )  where  S  =  average  stay  time  for  the 

type  unit  in  the  scenario. 

B-5.  ASSUMPTIONS  MADE  WITHIN  THIS  PAPER 

a.  For  a  unit  detection,  in  a  2-hour  period,  which  is  cued  by  a  detection  in  the  preceding  2-hour 
period,  Pr(detect  during  this  cycle  |  detection  previous  cycle]  =  y  •  POTA,  1  <  7  <  1/POTA.  We 
assume  that  y  has  a  Uniform  (1,  1/POTA)  distribution. 

b.  The  number  of  (whole)  2-hour  periods  after  initial  acquisition,  during  which  a  cued 
acquisition  is  retained,  is  distributed  as  Geometric  (1  -  y  •  POTA) 

c.  The  average  retention  time,  E[T|^]  ,  has  a  maximum  value  of  36  hours  in  NUFAM  because 
any  Tj^  exceeding  36  hours  is  equivalent  to  an  arbitrarily  long  time  in  a  NUFAM  simulation. 
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&-6.  ASSUMPTIONS  MADE  IN  NUFAM  lU  RELATING  TO  THE  ENGAGEMENT  OF  A  UNIT 
BY  NUCLEAR  WEAPONS 

a.  In  order  for  a  unit  to  be  available  for  engagement  (subject  to  fire)  by  nuclear  weapons,  it 

must  be  on  the  acquisition  list  at  the  end  of  the  fire  planning  time  (time  tp)  and  it  must  retained  as 
a  target  at  least  until  the  detonation  time  T,^.  We  define  the  probability  P[  unit  is  acquired  and 
retained  until  detonation]  as  Thus 

Pa»ai7~  ^  Tj  — tp  |  acquired  at  tp  ]  •  P[  acquired  at  tp  ]. 

Po»a«/~  Iq  >  *<l~^  I  Tj— tp  =  t  ]  P[  T,J  — tp  =  t  ]  dt  •  Pacq- 

b.  The  probability  that  a  unit  can  be  fired  at  by  a  nuclear  weapon,  denoted  as  is  the  joint 

probability  that  ::  unit  is  on  the  acquisition  list,  is  not  dropped  from  the  acquisition  list  prior  to 
detonation,  and  is  stationary  at  the  time  of  last  observation.  This  probability  calculation  implicitly 
assumes  that  a  weapon  of  the  appropriate  yield  is  avmlable  and  will  be  fired  at  the  target  if  it  meets 
these  criteria. 

Pfirtd~  P«««  ■  Iq  ^  ^d~*'P  I  t  ]  P[  Tj  — tp  =  t  ]  dt  •  p,^^^ 

where  Pi^f  =  P[  unit  was  stationary  at  the  time  it  was  last  observed  ]. 

c.  The  probability  that  a  unit  is  hit,  denoted  as  ,  is  the  joint  probability  that,  a  unit  is  on 
the  acquisition  list,  is  not  dropped  from  the  acquisition  list  prior  to  detonation,  is  stationary  at  the 
time  of  last  observation,  and  does  not  move  from  the  last  computed  OGZ.  This  calculation 
implicitly  assumes  that  a  weapon  is  available  for  firing  at  the  unit,  the  weapon  is  fired,  and  both  the 
weapon  system  and  warhead  function  correctly.  Thus 

PAiI  =  P«»oi7  ■  P»® 

Pfcfi  =  p«®«  ■  Iq  I  Tj-tp  =  t  ]  P[  T,,-tp  =  t  ]  dt  •  p,^,,  , 

PAi<  ~  Pfircd  '  g  .  vj  ’ 

where  is  computed  as  given  in  the  previous  section  and  SJ^^  and  are  the  values  entered 

into  the  NUFAM  III  data  base  for  the  unit  stay  time  and  move  time  respectively. 
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Section  ffl.  TARGETS  ENGAGED  WHILE  STATIONARY 
THAT  FACE  NO  TRACKING  CAPABILITY 


B-7.  SCENARIO.  The  type  unit  is  any  type  other  than  those  identifled  as  an  “engaged  while 
moving”  type.  We  assume  that  the  target  is  sufficiently  deep  relative  to  the  PLOT  such  that  there  is 
no  capability  to  continuously  (or  intermittently  with  a  short  duty  cycle)  monitor  or  track  the  target. 
The  target  is  acquired  while  stationary  or  is  acquired  in  movement  toward  an  identified  stationary 
location  (e.g.,  assembly  area),  where  d  represents  the  time  remaining  from  the  time  of  acquisition 
until  it  stops  at  this  stationary  location.  We  only  successfully  engage  the  unit  if  it  remains  in  place 
at  the  point  at  which  they  were  acquired.  We  assume  a  small  d  such  that  d  <  |  M. 


B-8.  ZONES.  NUFAM  III  zones  4  and  5. 


B-9.  STAY  TIME  AND  MOVE  TIME. 

S^  =  S  ,  the  scenario  average  time  the  type  unit  will  remain  stationary. 

M^=  M  ,  the  scenario  average  time  that  a  unit  will  be  moving  from  stationary  location  to 
stationary  location. 


B-10.  PROBABILITY  THAT  A  UNIT  WAS  STATIONARY  AT  THE  TIME  OF  LAST 
OBSERVATION.  This  probability  is  equal  to  the  probability  that  the  unit  was  stationary  at  the 
time  of  acquisition  or  that  it  was  moving  but  within  d  minutes  of  stopping  at  a  definable  stationary 
location.  Thus 


PROB.STAY.LOBS  =  +  S  •  -  ^  ^ 

S  +  M  S  +  M 


S  +  M 


d 

M 


B-11.  RETENTION  TIME.  The  retention  time  is  an  estimate  of  the  length  of  time  that  the  target 
unit  will  remain  at  the  location  at  which  it  was  acquired. 


R  =  ys 


+  s 


d  _ 
M 


■  14  ,  d  ■ 
.57  M.  • 


B-12.  ALTERNATE  REPRESENTATION.  Another  representation  does  not  use  the  probability 
that  the  unit  was  stationary  at  the  time  of  last  observation  (the  same  as  setting  PROB.STAY.LOBS 
equal  to  1  ).  It  is  not  quite  as  accurate  as  the  preceding,  but  does  not  require  the  code  associated 
with  a  check  of  PROB.STAY.LOBS.  As  a  result,  we  need  to  calculate  the  unconditional  retention 
time  expectation  E[  T^  ]  rather  than  the  conditional  expectation  E[  T^  |  stay  ].  The  formula  for 
E[  T«]i8: 

E[Tfl]  =  E[  Y^  I  stationary  at  time  of  acquisition  ]  •  P(  stationary  ] 

+  E[  S  I  arriving  at  a  stationary  location  at  time  of  acquisition  while  moving  ] 

•  P[  arriving  at  a  stationary  location  at  time  of  acquisition  |  moving  ]  •  P[  moving  ] 

Th„il  =  E(Y,].^  +  ElS|.«.  ']■ 

This  formula  should  only  be  used  if  PROB.STAY.LOBS  =  1. 
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Section  IV.  TARGETS  ENGAGED  WHILE  MOVING 
THAT  FACE  NO  TRACKING  CAPABILITY 


B-13.  SCENARIO.  The  type  unit  is  identified  as  a  “engaged  while  moving”  type.  We  assume  no 
capability  to  continuously  monitor  or  track  the  location  of  a  target  unit  acquired  as  part  of  a 
column  of  units  moving  along  a  specified  route.  The  nuclear  weapon  is  targeted  at  the  location  of 
the  unit  when  acquired  —  because  the  unit  is  one  of  many  moving  along  that  route,  it  is  anticipated 
that  a  similar  unit  will  be  within  the  radius  of  effects  of  the  weapon  when  it  detonates. 


B-14.  ZONES.  NUFAM  III  zones  4  and  5. 
B-15.  STAY  TIME  AND  MOVE  TIME. 


length  of  each  unit  in  column  (km) 
speed  of  column  (km/hr) 


length  of  gaps  between  units  in  the  column  (km) 
speed  of  column  (km/hr) 


NOTE:  Since  Mjy^  must  be  the  same  for  all  zones  (only  one  entry  is  allowed  for  each  type  unit),  Sj^^ 
for  each  type  unit  will  normally  be  the  same  for  all  zones  as  well. 


B-16.  PROBABILITY  THAT  A  UNIT  WAS  STATIONARY  AT  THE  TIME  OF  LAST 
OBSERVATION.  This  probability  is  not  applicable  in  this  situation;  we  ensure  that  NUFAM  III 
does  not  consider  it  by  setting  the  input  value  to  1. 

PROB.STAY.LOBS  =  1 


B-I7.  RETENTION  TIME.  The  retention  time  is  an  estimate  of  the  length  of  time  that  other  target 
units  will  be  moving  through  the  location  at  which  the  unit  was  2u:quired. 

_  (  Column  length  in  km  )  /  2 
speed  of  column  (km/hr) 
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Section  V.  TARGETS  FACING  A  NEARLY  CONTINUOUS  TRACKING  CAPABILITY 

(WITH  SCHEDULED  FIRE) 


E-18.  SCENARIO.  The  unit  is  located  relatively  close  to  the  FLOT.  We  assume  that  we  can 
continuously  (or  intermittently  with  a  short  duty  cycle)  monitor  or  track  the  location  of  an  acquired 
target.  The  desired  ground  zero  (DGZ)  of  the  weapon  is  shifted  along  with  the  unit  until  the  final 
adjustment  is  made.  Let  ti  denote  the  amount  of  time  between  the  final  observation  (time  of  final 
fire  adjustment)  and  the  detonation.  We  assume  a  small  q  <  min{  |  S,  |  ^  where  S  and  M  are 
defined  below.  From  another  perspective,  the  unit  (whether  stationary  or  moving)  remains 
stationary  wHh  respect  io  the  DGZ  during  this  entire  tracking  period.  The  time  of  detonation,  Tj,  is 
scheduled  (set)  during  the  time  of  fire  planning. 


B-19.  ZONES.  All  units  in  NUFAM  III  zones  I,  2,  and  3  engaged  with  scheduled  fire. 

B-20.  STAY  TIME  AND  MOVE  TIME. 

Let  S  =  the  scenario  average  time  the  type  unit  will  remain  stationary. 

Let  M  =  the  scenario  average  time  that  the  type  unit  will  be  moving  from  stationary  location 
to  stationary  location. 

Let  S^  and  M j^f  denote  the  stay  and  move  times  actually  used  as  NUFAM  III  input. 

17  denotes  the  amount  of  time  between  the  final  observation  (  time  of  final  fire  adjustment  )  and 
the  detonation. 

=  »? 

Sa,  =  S  -  Tj 


B-21.  PROBABILITY  THAT  A  UNIT  WAS  STATIONARY  AT  THE  TIME  OF  LAST 
OBSERVATION. 


PROB.STAY.LOBS  =  p,.,,  . 


B-22.  RJSTEN'llON  TIME.  The  retention  time  is  an  estimate  of  how  long  the  target  can  be  treu:ked 
using  the  same  sensors  that  initially  detected  the  unit.  Let  7  be  a  factor  that  indicates  how  much 
better  the  sensors  ar*  at  detecting  a  target  previously  detected  and  identified  than  detecting  a  new, 
previously  undetected  target  unit,  where  1  <  7  <  1/POTA,  and  all  values  of  7  between  1  and 
1/POTA  are  considered  equally  likely. 


POTA  <  0.97: 


R 


POTA  -  0.97 
(  1  -  POTA 


In 


[■ 


1  -  0.97 
1  -  POTA_ 


+ 


1.08 

1  -  POTA 


(in  hrs). 


POTA  >  0.97:  R  =  36  hrs. 
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Section  VI.  TARGETS  FACING  A  NEARLY  CONTINUOUS  TRACKING  CAPABILITY 

(WITHOUT  SCHEDULED  FIRE) 


B-23.  SCENARIO.  The  unit  is  located  relatively  close  to  the  FLOT.  We  assume  that  we  can 
continuously  (or  intermittently  with  a  short  duty  cycle)  monitor  or  track  the  location  of  an  acquired 
target.  The  desired  ground  zero  (DGZ)  of  the  weapon  is  shifted  along  with  the  unit  until  the 
beginning  of  the  fire  execution  cycle,  if  the  unit  is  stationary  at  the  time  that  the  fire  execution  cycle 
begins,  or  until  the  first  time  that  the  unit  stops  during  the  fire  execution  cycle,  if  it  was  moving  at 
the  time  of  the  beginning  of  the  fire  execution  cycle  and  it  does  stop  during  the  fire  execution  cycle. 
At  that  time,  the  time  of  the  detonation,  Tj,  is  scheduled.  Let  r;  denote  the  amount  of  time  between 
the  final  observation  (time  of  final  fire  adjustment)  and  the  time  of  detonation.  From  another 
perspective,  the  unit  (whether  stationary  or  moving)  remains  stationary  with  respect  to  the  DGZ 
during  this  entire  tracking  period. 


B-24.  ZONES.  All  units  in  NUFAM  III  zones  1,  2,  and  3  engaged  without  scheduled  fire. 


B-25.  STAY  TIME  AND  MOVE  TIME. 

Let  S  =  the  scenario  average  time  the  typ>e  unit  will  remain  stationary. 

Let  M  =  the  scenario  average  time  that  a  unit  will  be  moving  from  stationary  location  to 
stationary  location. 

Let  and  denote  the  stay  and  move  times  actually  used  as  NUFAM  III  input. 

T)  denotes  the  amount  of  time  between  the  final  observation  (  time  of  final  fire  adjustment  )  and 
the  detonation. 

L  denotes  the  length  of  the  fire  execution  cycle. 


SAf  =  P[Y5>i?]-  -^  +  P(S>q].P(Y,„<L-q]. 

Mjv,  =  (l-P(Y5>„]).  _l_  +  (i-P[S>q]).p[Y„<L-q]. 
where 

P[Y5  >  ,7]  =  1  -  0<  y<|s, 

0  <  u  <  I  S, 

|5<u<  ^s. 

0  <  L  — 7  <  5  M, 

2  M  <  L-7  <  I  M. 


P[YAf  <L-7l=  V 


M 

1 

M 


2(L-7)  - 


3  (L-7) 
4  M 


P[  S  >  7  ]  =  1  , 

—  9  _  3u 

E  ’ 
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B-26.  PROBABILITY  THAT  A  UNIT  WAS  STATIONARY  AT  THE  TIME  OF  LAST 
OBSERVATION. 


PROB.STAY.LOBS  =  --^5  +  P[  Ym  <  L-ij]  • 

S  +  M  S  +  M 

B-27.  RETENTION  TIME.  The  retention  tjme  is  an  estimate  of  how  long  the  target  can  be  tracked 
using  the  same  sensors  that  initially  detected  the  unit.  Let  7  be  a  factor  that  indicates  how  much 
better  the  sensors  are  at  detecting  a  target  previously  detected  and  identified  than  detecting  a  new, 
previously  undetected  target  unit,  where  1  <  7  <  1/POTA,  and  all  values  of  7  between  1  and 
1/POTA  are  considered  equally  likely. 


POTA  <  0.97: 

j,  _  POTA  -  0.97 
(  1  -  POTA  )= 


POTA  >  0.97:  R  =  36  hrs. 


In  r  1  1  +  _ _  (in  hrs) 

L  1  -  POTaJ  ^  1  -  POTA 


NOTE:  When  NUFAM  III  determines  if  a  target  unit  will  be  retained  at  least  until  the  time  of 
detonation,  the  P[  unit  retained  until  time  T ^  |  unit  on  acquisition  list  at  time  T j  ]  should  equal  1 
for  units  in  this  category.  The  current  code  does  not  distinguish  between  this  type  of  unit  and  the 
scheduled  fire  type  of  unit,  so  the  probability  P(  unit  retained  until  time  Tj  |  unit  on  acquisition  list 
at  time  Ty  ]  =  P[  Yr  >  nl  unit  on  acquisition  list  at  time  Ty  ].  However,  the  error  introduced  by 
this  calculation  should  be  very  small  since  R  tends  to  be  much  greater  than  t). 

B-28.  ALTERNATE  REPRESENTATION.  Another  representation  does  not  use  the  probability  that 
the  unit  was  stationary  at  the  time  of  last  observation  (the  same  as  setting  PROB.STAY.LOBS 
equal  to  1  ).  It  is  not  quite  as  accurate  as  the  preceding,  but  does  not  require  the  code  associated 
with  a  check  of  PROB.STAY.LOBS.  As  a  result,  is  calculated  in  the  same  way  as  given  above 
but  Mjjy  =  1  —  S^: 


Sjv,  =  P[Y5  >  »,]  .  +  P[S  >  r,]  ■  P(Y„,  <  L-r,]  • 


M  _ 
S  +  M’ 


M  w  =  1  -  Sy 


■ 
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APPENDIX  C 

PROBABILITY  DISTRIBUTIONS 


C-1.  GEOMETRIC  DISTRIBUTION 
Suppose  that  X  ~  Geoinetric(  p  ).  Then 

P[X  =  r]  =  p(  1-p)*,  z=0,  1,  ... 

The  expectation  E[  X  ]  = 

C-2.  STAY  TIME  DISTRIBUTION 
Uniform  (  j  S,  |  S  ),  where  S  is  the  Average  Stay  Time. 
P[S>u]  =  F5(u)=  1, 

is-  is  “  2S’ 


0  <  u  < 

|S  <u<  Is. 


C-3.  RESIDUAL  STAY  TIME  DISTRIBUTION 

Let  Y 5  denote  the  residual  life  of  the  stay  time.  Then 

P[V,>»1  =  1-  EPl  lo '  -  I 


=  1 

y 

s  ’ 

0  <  y  < 

_  4 
“  3 

1 

S 

2^-  '4] 

4  S  J 

|S  <  y 

-  s 

3 


4  s 
3 


C-4.  EXPECTATION  OF  THE  RESIDUAL  STAY  TIME 

EiYji  =  1“  y  •  fy^fy)  '■y  =  ir  ^ 

1  S  3  S 

=  E^  1  y  •  ^  +  gfe]  1  -y  ■  1" 


_  14  c 

-  5=s. 


§5 
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C-5.  MOVE  TIME  DISTRIBUTION 
Uniform  (  |  M,  5  M  ),  where  M  is  the  Average  Stay  Time, 


Fm  ('»)  = 


1«- 


4  O 


M  -  Tm 


_  9  3u 
2M 


i 


<  u 


<  ^M. 


C-6.  RESIDUAL  MOVE  TIME  DISTRIBUTION  AND  EXPECTATION 

Let  Y nf  denote  the  residual  life  of  the  move  time.  Since  S  and  M  have  the  same  distribution,  so 
do  Y5  and  Yjy^;  simply  replace  S  with  M  in  the  formula  in  paragraph  C-3. 


C-7,  COMPUTING  P[  Y^  >  Tj  -  T.c,  ] 

The  time  Tae^  is  not  known,  but  the  difference  (  tp  —  Tac,  )  is  distributed  as  the  age  of  the 
distribution  of  the  time  of  acquisition,  Y^^  (the  age  and  the  residual  life  are  identically  distributed), 
independent  of  the  distribution  of  Tj  — tp.  Tj— tp  can  have  any  nonnegative  distribution;  in  this 
case,  a  uniform(  0,  L  )  distribution  is  assumed,  where  L  is  the  length  of  the  fire  execution  cycle. 
Thus  ^ 

P(  ^5  >  Tj  —  Taej  1  =  I  |p-o|,  g  — tp)  +  (tp  — Taej)  |  Tj  — tp  =  t;  (tp  — T«cg)  =  y] 

■  F(Tj— tp=:t]  P[tp— Tocg  =  y]  dt  dy 

where  D*  =  min{  L,  j  S  } 


PIY5  >  t+y]  =  1  - 


t+y 

s  ’ 


0  <  t  +  y  <58, 


P[T,-tp  =  t] 
P[tp— Tacj=y] 


_  4 
“  3 


_  1 
~  L 


S 


2(t+y)  - 


3  (t+y)^ 

4  S 

0  <  t  <  L. 


,  g  S  <  t+y  <  I  S. 


-  >  y  1  _  I  0  <  V  <  R 

“  E[  Tfl  ]  ~  a  ’  ®  y 


NOTE:  Since  (Tj— tp)  and  (tp— Tac?)  are  both  distributed  as  unifoim  random  variables,  it  is  easier 
computationally  to  determine  the  distribution  of  W  =  (Tj  — tp)  +  (tp— T^c,)  and  calculate 


P[  Y5  >  Trf  -  Tac,  I  =  P(Y5>  W  (  VV  =  w  j  P(  W  =  w  ]  dw 
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